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!>• 1 ABSTRACT 

We present a general overview and the first results of the SUMO project (a SUrvey of Multiple 
populations in Globular Clusters). The objective of this survey is the study of multiple stellar 
populations in the largest sample of globular clusters homogeneously analysed to date. To 
CO , this aim we obtained high signal-to-noise (S/N>50) photometry for main sequence stars with 

■ mass down to ~0.5 Alodot in a large sample of clusters using both archival and proprietary U, 

CO ' B, V, and / data from ground-based telescopes. 

In this paper, we focus on the occurrence of multiple stellar populations in twenty three clus- 
ters. We have defined a new photometric index cu,b,i - (U — B) - (B - I), that turns out to be 
very effective for identifying multiple sequences along the red giant branch (RGB). We found 
that in the V"-ctj,b,i diagram all clusters presented in this paper show broadened or multimodal 
RGBs, with the presence of two or more components. We found a direct connection with the 
chemical properties of different sequences, that display different abundances of light elements 
(O, Na, C, N, and Al). The cu,b,i index is also a powerful tool to identify distinct sequences of 
stars along the horizontal branch and, for the first time in the case of NGC 104 (47 Tuc), along 
the asymptotic giant branch. Our results demonstrate that i) the presence of more than two 
stellar populations is a common feature amongst globular clusters, as already highlighted in 
previous work; ii) multiple sequences with different chemical contents can be easily identified 
by using standard Johnson photometry obtained with ground-based facilities; Hi) in the study 
of GC multiple stellar populations the ctj.b.i index is alternative to spectroscopy, and has the 
advantage of larger statistics. 
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Based on observations made with MPG/ESO 2.2m telescope in La Silla, 
Chile, proposal 088.A-9012. 
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1 INTRODUCTION 

There is an increasing body of evidence suggesting that both Galac- 
tic and Magellanic Cloud globular clusters host multiple popula- 
tions, a byproduct of the cluster internal chemical evolution. The 
majority of GCs appear to have at least two approximately coeval 
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sub-populations, or rather two generations of stars. The recognition 
of this fact came with the increasing observational evidence for 
what was denoted as 'star-to-star abundance variations' or 'abun- 
dance anomalies', that manifest themselves as variations of nitro- 
gen, carbon, oxygen and sodium - in some clusters also magnesium 
and aluminum - abundances amongst stars in individual clusters. 
These abundanc e differences are found in stars of all evolutionary 
stage s (see, e.g. iKraft et alJl992l ; ICannon et aljl9 98; Gratt on et al] 
1200 lh . 

This phenomenon ha s been known for a few decades (see, e.i 
IOsbornlll97ll; iKrafdl 19781 : ICohenlll978l : iNorris & Free"manlll979 
but was limited to a small number of stars in each cluster. The strik- 
ing, additional evidence that GCs host multiple populations came 
from both high-quality photometry for thousands of cluster stars, as 
well as from systematic spectroscopic investigation thanks to mod- 
ern multi-object spectroscopy. As it turned out, multiple cluster se- 
quences in the colour-magnitude diagrams (CMDs) show up in all 
the clusters studied with appropriate filters. The occurrence of mul- 
tiple populations can be related to th e evidence of m a in sequence 
(MS) splittings like in a) Centauri ( I Anderson 19971: Bedin et alj 



|2004) . NGC 2808 dD'Antona et alj|2005l ; IPiotto et alj|2007l) . or 
NGC 6752 dMilone etal.ll201Cl) . but it is also associated to the 



photometric detection of multimodal RGBs (e.f 


r. Grundahl et alj 


ll999l;lMarino et alj2008l;lYone & Grundahll2008 


;lLee et al.l2009). 



and in some cases of a sub-giant branch (SGB) splitt ings (e.g. 
iMilone et al]2008l : lMarino et alj|2009l : |piotto et alj|2012h . 

The most recent photometric analyses have shown that the 
CMDs of GCs consist of intertwined sequences of the two or more 
populations, whose separate identities can be followed continu- 
ously from the MS up to the red giant branch (RGB), and t hence to 
the horizontal branch (HB, NGC 104, IMilone et alj2012el) . 

Photometric analyses have clearly shown that the use of pho- 
tometric filters sampling the ultra-violet (UV) portion of the stel- 
lar spectra, is a very powerful tool for detec ting and tracing mul- 
tiple populations in the CMDs of GCs (e.g. iGrundahl et all 19991 : 
iMarino et alj|2008h . 

Computations of suitable model atmospheres ( Sbordo ne et alj 

1201 lh accounting for the peculiar chemical patterns of second gen- 
eration stars have explained the pivotal role played by UV photo- 
metric filters in studying sub-populations in GCs. The UV portion 
of stellar spectra is hugely affected by the change in the abundance 
of CN and NH molecules, due to the light-element anticorrelations 
observed in the stars belonging to the second stellar generation; 
using suitable combinations of photometric filters, including UV 
filters, allows one to emphasize the difference in the radiative flux 
between primordial cluster stars and those enriched by the high- 
temperature H-burning products. 

Obviously, (high-resolution) spectroscopy allows a more de- 
tailed analysis of the specific chemical patterns associated to sec- 
ond generation stars; but for the majority of GCs this kind of analy- 
sis is limited to the brightest RGB stars, and in any case to a limited 
number of objects. 

Tracing the multiple population phenomenon in a large frac- 
tion of cluster stars (if not in the whole cluster) is mandatory to es- 
tablish the actual fraction of second generation stars in individual 
GCs. In addition, the comparison betw een photometry and th eoret- 
ical stellar evolutionary models (e. g. IPietrinferni et al.ll2009h and 
colour - r cff transformations, that account fo r the specific chemica l 
patterns of the distinct sub-populations (see ISbordoneetalj201lh . 
allows to investigate how stellar structural and evolutionary prop- 
erties are affected by the chemical abundance distribution of each 
distinct sub-population. 



In this work we present an overview of the SUMO project, - 
a SUrvey of Multiple populations -, together with its first results. 
The paper is organized as follows: Sj2] introduces the main objec- 
tives of the survey, while the data and the derived CMDs are pre- 
sented in ^3]and {j4] ^introduces a new photometric index, based 
on the combination of three bands, that is effective to highlight dif- 
ferent sequences along the RGB. The connection with the chemical 
content of sub-populations identified in this way is discussed in ijfj] 
Our conclusions in 0close the paper. 



2 THE SUMO PROJECT 

The basic idea of the SUMO project is to create a database of ho- 
mogeneous, wide-field, multiwavelength, photometric data of GCs. 
The main aim is to study the multiple population phenomenon, and 
in particular the primary objectives of the survey are: 

• investigate the occurrence of multiple populations in a large sam- 
ple of GCs. Is this a property common to all GCs? 

• Identify multiple populations along different evolutionary phases, 
from the MS to the RGB and HB; 

• verify the connection between photometric and spectroscopic 
properties of different stellar populations; 

• study the radial gradients over a large fraction of the body of each 
cluster. 

U,B,V,I photometry for a large number of clusters will also 
be used to study a number of different topics, to be presented in 
forthcoming papers: 

i) the effect of diffe rential redden i ng, if present, following the pro- 
cedure described in IMilone et alj d2009l see also Sect. 14. 1) ; 

ii) the population of blue stragglers, and their radial gradients; 

Hi) the binary fraction, which will complement similar works on the 
same GCs based on Hubble Spa ce Telescope (HST) data and lim- 
ited t o the cluster centres (e. g. ISollima et alj 120071 : IMilone et alj 
l2008h . 

iv) the properties of the RGB bump over a large metallicity range; 
vj a self-consistent study of the HB morphology; 

vi) identification of spectroscopic targets for future follow-ups; 

vii) new theoretical models will have to be developed, with ad-hoc 
chemical mixtures an d the proper model at mospheres, following 
the pioneering work of lSbordone et aljj201lh . 

This survey has been designed to complement the HST data 
available in the literature, that have been successfully used to dis- 
entangle multiple sequences in the CMP of the inner, most crowded 
regions of a large numbe r of GCs (Piotto et al. 2007; Milon e~et alj 
l2008l: iBelliniet all l20ld : IMilone et alj|2012el : IPiotto et alj|2012h . 
The large field-of-view of ground based data is mandatory to extend 
this result over the main body of each cluster, tracing the presence 
and the relative contribution of distinct components. This is crucial 
to shed light on the formation mec hanism of the multip le popu- 
lations. In fact, theoretical models dP'Ercole et al. predict 
that a large fraction (up to 90% or more) of first-generation stars 
should be lost early on during the GC evolution, because of the 
expansion and stripping of the cluster outer layers resulting from 
the mass loss consequent to massive SNe explosions. The fact that 
first-generation stars were much more numerous at the time of the 
cluster formation can also account for the polluting material needed 
to explain the following stellar generation. Signatures of these dy- 
namical evolution are expected to be observable in the present-day 
stellar radial distribution, at least in clusters with long relaxation 
times; a systematic study of the radial gradients would therefore 
provide solid constraints to understand the nature of the polluters 
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Figure 1. Left-Reddening map for NGC 6121. Center and Right- V vs B - I CMD before and after applying the correction for differential reddening. 



and the processes that have caused the formation of successive stel- 
lar generations. 

Finally, follow ing previous works llMarino et all [2008; 
iMilone et al1l2012eh . we emphasize the use of the {/-band as an 
efficient diagnostic to identify multiple stellar populations. For this 
reason, from the observational point of view (see i}3}, we devoted a 
significant fraction of our observing runs collecting data in this fil- 
ter. The colour combination presented in this paper (see §.[5]( turned 
out to be very effective to show the presence of multiple popula- 
tions along the RGB, correlating with the chemical properties of 
stars. In the following section, we summarize the properties of the 
target clusters. 



2.1 Public products 

One of the objectives of this project is to disseminate the results 
to the greatest extent possible, releasing to the general commu- 
nity a large number of products. With this aim, we set a web 
page, http://www.iac.es/project/sumo which will be regularly up- 
dated with: plots (CMDs, colour-colour), reddening maps, fiducial 
lines. 



2.2 Properties of the target clusters 

The cluster sample presented in this paper includes both objects 
where different generations of stars have been already discovered 
and widely studied, and objects with no evidence of multiple stellar 
populations. In this Section we summarize the observational sce- 
nario. The main properties will be used in the following to interpret 
the multimodal RGBs of Fig. [6TTT1 

The CMD of 47 Tuc-NGC 104 consists of intertwined se- 
quences, whose separate identities can be followed continu- 
ously from the MS up to the RGB, and thence to the HB 
(Milo neetai]|2012eh . The observed colours are consistent with 
a pa ir of populations with different content of light elements 
(e.g. iNorris & Freeman! 1 1979t ICannon et all 1 19981: lHarbeck et all 
2003r" 



(Gratt on et at] 1 19861: IBrown et al.1 Il99fi iBrow n & Wallerstein 



19921 : lDrakeet al.il 19921; ISmith et al.1 120051: 



.. IVillanova & Geislet 

l201ll) with a bimodal CN distributi on jNorrisll98lh. an O-Na anti 



Pancino et al. 2010h. and small helium differences ( e.g. 



di Criscienzo et alj201 ll : lNataf et al.l2"oTT| ; lMilone et alj2012dh . A 
third population is vi sible only along the SGB, and includes less 
than 10% of the stars dAnderson et al.ll2009l : |piotto et alj|2012h . 

NGC 6121 is known to be a metal-intermediate GC, with 
[Fe/H]a -1.1 and [or/Fe]«0.4. Spectroscopy of bright RGB 
stars revealed a significant spread in light-element abundances 



corre lation and a Na-Al correlation divans et al.ll999tlMarino et alj 
2008). The analysis of high-resolution VLT/UVES spectra for 
about one hundred RGB stars reveale d a bimodality in the stellar 
distribution in the Na versus O plane (Mari no et al . 2008). The two 
groups of Na-rich/O-poor, and Na-poor/O-rich stars populate two 
different sequences along the RGB in the U versus U - B CMD. 
Na-rich stars define a sequence on the red side of the RGB, while 
Na-poor stars populate a bluer, more spread sequence. Interestingly, 
NGC 6121 hosts also a bimodal HB, which is well populated both 
on the red and the blue side of the RR-Lyrae gap. HB stars exhibit a 
bimodal Na and O distribution similar to what found for RGB stars, 
with red-HB stars having solar N a and blue-HB sta rs being oxygen 
depleted and sodium enhanc ed dMarino et alj|201 H see also lNorrisI 
ll98ll ; ISmith & Norrisl 19931 

Among clusters with multiple stellar populations, NGC 2808 
is certainly one of the most intriguin g objects. Its CMD shows 
three distinct MSs dPiotto et af] |2007|) . with the middle and the 
blue MS being highly enhanced in helium compared to the red 
MS - AY~0.07 and AY-0.13, respect ively dD'Antona etai1l2005l; 
IPiotto et alj|2007l : IMilone et al] |20 1 2d ) . It also hosts a multimodal 
HB g reatly extended blue -ward j [sosin et al.l [l997l ; iBedin et alj 
2000), and a spread RGB iLee et al.l d2009h . Furthermore, spectro- 
scopic studies of RGB, HB, and bright MS stars have revealed 
significant star-to-star differences in th e light-element abundances 
and an extended Na-0 anti correlation (Norris 19 8lt ICarretta et ail 
l2006l:lBragaglia et alfeoid) . 

In case of NGC 6205 large ranges of light-element abun- 
dances have been found along all evolutionary phases from the 
RGB tip down to the MS (e. g. IShetrondfl996l : iKraft et alJI 1991 
iBrilev et alj|2002h . NGC 6205 exhibits a very extend ed Na-O and 
Al-Mg anticorrelation, and a C-N correlation (see e. g. ICohen et alj 
120021 ; Isneden et al.ll2004l ; ISmolinski et alj|201ll and references 
therein). The RGB is widely spread in c olour in the y versus 
c v CMD, with at least three components ( Grundahl et al.1 1 19981 : 
lYong & Grundahll 120081 ; lLardo et al.ll201ll) . confirming that this 
cluster is made of stellar populations with different nitrogen and 
carbon content. 

Since the seventies, spectroscopic studies of NGC 6397 gi- 
ants have shown that the RGB stars exhibit a significant spread 
in C, N, Na, and O abundances (e.g. iBell et al" Hl979l ; lBrilevetalJ 
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Table 1. Summary of data for the 23 clusters included in this papers. We summarize the cluster name, coordinate, and the maximum number of images, in 
each band, a given stars has been measured in. 
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Figure 2. (U - V) vs (B - I) colour-colour plane for the NGC6121 (left) and NGC 6205 fields (right). Following iBono et all feplOl) we use this diagram to 
identify foreground Galactic stars and background unresolved galaxies. The bona-fide cluster stars are selected inside the area marked by the overplotted line. 
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Figure 3. Selection of different CMDs of NGC6121 and NGC6205. The plot shows that our photometry reaches well below the MS turn off, with good 
accuracy in all four bands. AGB and HB stars are highlighted with different colours to show their behaviour in the CMDs. 



ll99d : |Pasauini et alj|2004l; ICarretta et al.ll2009l) . and a bimodality 
in sodium and oxygen jLind et al" 1 120091 1201 lb . Multi-band HST 
photometry revealed that the MS of NGC6397 splits into two 
compo nents made up of ~ 30% and ~70% of MS stars, respec- 
tively dMiloneetal]|2012bh , while a bimodal RGB was previously 



identified by iGrundahl et all d2002l) and iLind et al.l d201 lb . Mul- 
tiple stellar populations along the HB have been investigated by 
lLovisi et al] d2012T) . The small colour separation of the two MSs 
implies that the two stellar populations differ in helium by AY~0.01 
ddi Criscienzo et alfcOlfjlMilone etafll2012bh . 
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Figure 4. NGC 6121 and NGC 6205 are adopted to illustrate the properties of the (V, cu.b.i) pseudo-CMDs. The left and central panels show the entire diagram 
and a zoom in the RGB region, respectively. Using this diagram, we identified two and three subpopulations, respectively, highlighted in the right panels. 



Multi-band HST photometry reveals that the MS, SGB, 
and RGB of NG C 6752 splits into three main components 
l lMiloneetal.ll2013h . in close analogy with the three distinct seg - 
ments along its horizontal branch stars (Mom anv et ail20 02. 2004). 
These triple sequences are consistent with three stellar groups: 
a stellar population with a chemical composition similar to field 
halo stars (population a), a population (c) with enhanced sodium 
and nitrogen, depleted carbon and oxygen and enhanced helium 
abundance (AY -0.03), and a population (b) with an intermedi- 
ate (between population a and c) light-element composition, and 
slightly helium enhanced (AY ~0.01). These components contain 
~25% (population a), ~45% (po pulation b), and ~30% (population 
c) of the cluster s t ars (see e.g. lNorris|[ l981; Grundahl et al. 2002; 
Yong et alj 120051 : ICarretta et al.l 120071: lYong & Grundahll booa 



Carretta et al I2012L" for spectroscopic studies on stellar populations 



in NGC 6752). 

In NGC 7 078 CN-CH C-N and Na-O anticorrelations have 
been observed ( Sneden et al .ll997tlCohen et alJbOOSLlkavser et alj 
120081 : IPancino et alj|2010[) and also a spread in the slow neutron - 
capture elements dSneden et alj 1 19971 12000| : ISobeck et all l201ll) . 



The RGB is broad ened in the u—g colour from SLOAN photometry 
jLardoetalfeoill) . 

In NGC 288, NGC 362, NGC 3201, NGC 4590, NGC 5904, 
NGC 6218, NGC 6254, NGC 6712, NGC 6809, NGC 7099, and 
NGC 6934 multiple stellar populations ar e inferred by star- 



to-star light-element abundance variations dPancino et 
"^n i„„„ j : — n u„„„i m 



5 



2010; 



2008; 



[ Cohen et al.1 120051: ICarretta et alj 120091; iMartell et all 
iLee et alj2009l : lKravtsov et alJl2009l : ISrmth & Bellll 19861) . 
HST photo metry reve a led a multiple sequences along the SGB o f 
NGC 362 IPiotto et alj J2012T) . As for NGC 288 dRoh et alJl2oTTh . 
and NGC 3201 the RGB is multimodal or spread when observed 
in appropriate UV bands. In NGC 6366, NGC 6093, NGC 6541, 
and NGC 6981 there is no evidence of multiple stellar populations 
to date, due to lack of spectroscopic studies of light-element abun- 
dances. 
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Figure 5. The ridge lines of the three populations detected along the RGB of NGC 6205 are used to show the efficiency of the ctj.b.i index in identify multiple 
sequences. The top three panels show the ridge lines in three different CMDs. The bottom panels show the colour difference with respect of the green line, 
calculated at V = 14.82 mag and V= 16.63 mag for the two clusters. The (U - B)- (B - I) index maximizes the colour difference, thus best disentangling the 
different sequences. 



3 DATA ACQUISITION, REDUCTION, AND 
STANDARDIZATION 

To accomplish the goals listed previously, we started a long-term 
observing campaign using facilities in both the northern and south- 
ern hemispheres, namely the Isaac Newton Telescope (INT) op- 
erated by the ING in La Palma (Spain), and the 2.2m MPG/ESO 
telescope located in La Silla (Chile). Both telescopes are equipped 
with comparable large field of view cameras: the WFC - Wide Field 
Camera (34'x34\ pixel scale 0.333"/px) and the WFI - Wide Field 
Imager (34'x33\ 0.238"/px), respectively. Both telescopes have a 
set of Johnson filters, including the {/-band one. 

The new data were used to complem ent the extensi ve data base 
maintained by one of us (P.B. Stetson, see IStetsonl2000l) , and in par- 



ticular to add the images in the U -band for a large number of clus- 
ters. So far, the data have been collected in four different observing 
runs at the INT, between August 2011 and August 2012, and one 
observing run at the 2.2m in February 2012 (complemented by few 
images collected in service time). We emphasize that this is a long- 
term, still ongoing project: the data reduction of the new data has 
been completed for a subset of the targets only, and first results are 
presented in the following sections. 

The pre-reduction was performed with standard IRAF proce- 
dures to correct for bias, flat-field and fringing in the case of the 
/-band images. The photo metric reduct i on wa s performed using 
DAOPHOT/ALLFRAME lStetsonlll987lll994l) simultaneously on 
all the data available for each cluster. Individual PSFs were mod- 
eled for each image, sampling the whole area of the chip with a 
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Figure 6. M v v.vcu,b,i C MDs for NGC 6366, NGC 104, NGC2808, and NGC 6121. The split of the RGB is evident all the clusters. We adopted the distance 
and mean reddening from [Harris to shift the diagrams to absolute magnitudes. Apparent magnitude and colour are labeled in the right and top axes, 

respectively. The clusters are sorted in order of increasing metallicity. 



large number stars (tens to hundreds), the exact number depend- 
ing on the stellar density of each field. After the ALLFRAME run, 
resulting ma gnitudes were corrected aperture using DAOGROW 
JStetsod[i99Ch . using a cons t ant co rrection. The calibration was 
performed following IStetsonl d2000b . So far, artificial stars tests 
have not been performed. Table [T] presents a summary of the data 
relative to the clusters included in the present paper. For e ach tar- 
get, w e indicate R.A. and Dec. (from the compilation by iHarrisI 
201G), and the maximum number of images a given stars has been 
measured in. 



4 PHOTOMETRIC SELECTIONS AND 
COLOUR-MAGNITUDE DIAGRAMS 



In this section we describe the procedure followed to obtain the 
cleanest possible catalogue for each cluster. Stars have been se- 
lected according to the photometric error, x 2 , and sharpness pa- 
rameters provided by ALLFRAME. In some cases we removed the 
stars located in the innermost, most crowded regions. The exact se- 
lections were optimized for each cluster. 



In the following subsections, we describe the technique 
adopted to correct for differential reddening and decontaminate 
for the presence of foreground and background Galactic stars and 
background unresolved galaxies. 
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Figure 7. Same as Fig|6] but for NGC 362, NGC 5904, NGC 288, and NGC 6218. 



4.1 Differential reddening 

We were careful about correcting our photometry fo r differential 
redden ing effects by using the procedure described in Milo ne et al.l 
(2009). Briefly, we first defined the ridge line for the cluster MS 
and RGB. For each star we identified the closest neighbours (typ- 
ically 50), and for each of them we estimated the amount of shift 
to the red or the blue side of the fiducial sequence. This systematic 
colour offset is indicative of the local differential reddening. We ap- 
plied this method using two CMDs: V vs (U - V) and / vs (B - I), 
to have two estimates using independent colour combinations. The 
final reddening map, presented in Fig.[T]for the case of NGC 6121 
(M4), was obtained as the mean of these two ones. The left panel 
shows the reddening map as a function of the XY coordinates, ac- 
cording to the colour code displayed at the top of the plot. In this 
case, we estimated that the extinction differences are of the order 
of AE(B - V) = ±0.05 mag. The central and right panels of Fig.Q] 



show the B vs (B - I) CMD of the same cluster before and after 
applying the correction for differential reddening. Clearly, all the 
features of the CMD become better defined, and the colour spread 
significantly decreases, especially in the case of the SGB, RGB and 
HB. Photometry corrected for differential reddening has been used 
in all the diagrams presented in the following. 

4.2 Removing contaminating sources 

To clean the CMDs from both foreground Galactic stars and back- 
ground unresolve d galaxies we followed a procedure similar to 
lBonoetalJ l l2010l) and based on the source position in the (B - 1) vs 
(U - V) colour-colour plane. Figure|2]shows the case of NGC 6205 
and NGC 6121. The sources inside the overplotted line are con- 
sidered bona-fide cluster stars, while those outside will be rejected 
from the following analysis. We stress that this procedure does not 
work equally well for every cluster. In fact, in case of more metal- 
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Figure 8. Same as Figg] but for NGC 6981, NGC 6934, NGC 6205, and NGC 6752. 



rich clusters, the RGB gets redder in the U—V colour, making more 
difficult to establish the separation from field stars. This is clearly 
seen in case of NGC 6121, where the non-negligible contamination 
from bulge stars cannot be completely removed, due to the super- 
position of the sequences in these diagrams. 

4.3 Colour-Magnitude Diagrams 

Figure [3] shows a number of different CMDs of the clusters 
NGC 6121 (left column) and NGC 6205 (right), adopted again as 
a test case, using different colour combinations. In particular, the 
top row presents the V vs (5 - /) diagrams, while the other panels 
display U vs (U - B), (U - V), and (U - 1), respectively. The figure 
discloses that our photometry reaches ~3 mag below the turn-off in 
the four bands. The sources rejected using the colour-colour plane 
are represented in light grey, while the bona-fine cluster stars are in 
black. We highlighted the AGB stars in orange, while HB stars are 



coloured according to the position relative to the " Momany jump" 
dMomanv et alj |2002|) and the "Grundahl jump" I Grund ahl et alj 
1999). In particular, blue, cyan, and green symbols represent the 
stars bluer than the Momany jump, intermediate, and redder than 
the Grundahl jump, respectively, while red dots mark the red HB 
stars. The selections were performed in the V vs (B - 1) plane. Note 
that in case of M 4, the red and the blue HB stars (red and green 
points, respectively) are swapped in the (17 - B) plane. Similarly, 
the bluest AGB stars (orange points) gets bluer than the reddest HB 
stars in the same plane. This behaviour is predicted by theoretical 
models, and is due to the sensitivity of U filter to the Balmer jump 
feature. 



11 



-2 

t — i — i — | — i — r - 

, H 

H . 

» * »^ 



U.B.I 

■1.9 -1.8 -1.7 
— | — i — i — i — i — | — i — i — i — r 



l — i — i — I — i — i — i — i — I — i — i — r 

.NGC6264 
[Fe/H]=-1.56 



- 14 o 



g _ • J $w- :if *v< ■ 




-2 



■1.2 



-1.1 



U.B.I 

-1 -0.9 -0.8 



-i — i — I — i — i — i — i — I — i — i — n — I — i — i — i — i — I — i — 

™ < NGd320i 
. ; • [Fe/H]=-1.59 



* o e 



- 14 




-i — I — i — i — i — i — i — i — n — i — I — rn — i — i — i — i — i — r 



Ho 
H 



_H 



NGC6681 
[Fe/H]=-t.62 



2 -a 



-2 



-0.9 -0.8 -0.7 
-1.9 -1.8 -1.7 



14 



H 

-« v 




T 



.T— i — i — I — i — r 



1 : Npd6093 ' -| i 



•=-> ° - 16 



U.BJ ^ U,B,I 

Figure 9. Same as Figg] but for NGC 6254, NGC 3201, NGC 6681, and NGC 6093. 



5 TIIK < | p. | INDEX AS A NEW DIAGNOSTIC FOR 
MULTIPLE STELLAR POPULATIONS 

5.1 Definition 

The selected catalogues were employed to look for the evidence 
of multiple stellar populations in the observed clusters. Recent 
studies have shown that CMDs of globulars consist of inter- 
twined sequences of two or more stellar populations, that can 
be t raced continuously f rom the MS up to the RGB and the 
HB jMiloneetalJl2012ij [cl). Previous stud ies of stellar popula- 
tions in GCs based on Stromgrem (e. g. Grunda hletal.ll 19991: 
lYong & Gru ndahl 200 ^), multi- wavelength HST and ground-based 
photom etry dMarino et alj|2008t iBellini et aLll2010l : iMilone et alj 
l2012al lcl) have demonstrated that appropriate combinations of 
colours and magnitudes can be powerful techniques to isolate mul- 
tiple sequences in the CMD of GCs. 

The (U - B) colour is sensitive to light-element variations 



dMarino et all 120081 : ISbordone etail l201ll) . These authors have 
shown that, at least in some clusters, stars enhanced in N and Na 
and depleted in O and C, populate the red part of the RGB in the 
U versus (U - B) CMD, while Na/N-poor and O/C-rich stars oc- 
cupy the red RGB region. A visual inspection at these diagrams 
reveals that the RGB is much broader than expected from photo- 
metric colour errors only, that for RGB stars are typically smaller 
than 0.03 mag. 

The (B - /) colour is very efficient in dis entangling stellar 
populations with differen t helium abundance JPiotto etall 120071 ; 
Idi Criscienzo et all 1201 lh . Both theoretical arguments and obser- 
vations show that Na/N poor stars are also helium-poor while stars 
enhanced in sodium and nitrogen have also higher helium abun- 
dance. Moreover, helium-rich stars have bluer (B - I) but redder 
(V - B) colour than helium-poor ones. Driven by these results, we 
introduce here a new index cu,bj, defined as (U - B) - (B - 1), that 
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Figure 10. Same as Fig© but for NGC 6541, NGC 6809, NGC 6397, and NGC 4590. 



maximizes the separation among stars with different helium and 
light-elements content. 

In the following, we will use NGC 6205 and NGC 6121 to de- 
scribe the properties of the cu b.i index. Figure |4]presents the V vs 
c u,b,i diagram for both clusters: the left and central panels show 
the entire diagram and a zoom in the RGB region, respectively. In 
this plane, the pseudo-CMD appears reversed when compared to a 
typical CMD: the c V B1 index gets larger (i.e. less negative) for MS 
stars of increasing brightness, while the SGB bends towards more 
negative colours. Interestingly, in the case of NGC 6121 the red HB 
crosses the RGB and the blue HB extends to the larger cu b.i index, 
while in the case of NGC 6205 it is completely outside the range 
of the plot (cu.b.i > -0.5mag, and not shown for clarity. An homo- 
geneous investigation of the HB morphology will be presented in a 
forthcoming paper. 

Contrary to what shown in Fig. [3] the RGB shows a large 
colour spread, significantly larger that the photometric error, and 



multiple sequences can be recognised in both clusters. This sug- 
gests that multiple stellar populations are present and that the 
c u,b,i index is able to identify them, by splitting the RGB into 
multiple components. In particular, in the case of NGC 6121 Fig- 
ure [4] suggests that the RGB is bimodal in the magnitude interval 
~ 12.5 < V <~ 16.1. We have thus drawn by hand the grey lines of 
Fig. E] (top right panel) to separate two groups of stars in the mag- 
nitude bin where the RGB split is more evident (12.5 < V < 16.1). 
These stars are coloured green and magenta, and named RGBa 
and RGBb, respectively. Also in the case of NGC 6205 the RGB 
is broadened, with the presence of at least three components. Sim- 
ilarly to what done for the case of NGC 6121, in Fig. |4] (bottom 
right) we have selected three groups of stars along the RGB of 
NGC 6205 (RGBa,b,c), that have been plotted with green, magenta, 
and cyan colours. 

The selected populations help to understand why the c U B ,i is a 
good colour combination to split the different RGBs that otherwise 
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Figure 11. Same as Fig|6]but for NGC 7099 and NGC 7078. 




Figure 12. Sketch of the procedure to estimate the RGB width for NGC 6205. Panel (a). V-cu,b,i diagram zoomed around the RGB. The red-dashed line is 
the RGB fiducial line and is drawn by hand. Panel (b). Normalized V-Acu b,i diagram for RGB stars between the two grey lines of Panel (a). See text for 
details. Panel (c). Histogram of the colour distribution for the RGB stars displayed in panel (b). 



are photometrically degenerate. This is shown in Fig. [5] The top 
three panels show the ridge lines of the two sequences detected in 
NGC 6121, in three different CMDs (V vs (U - B), (B - V), (B - 1)). 
The same is shown for the three sequences of NGC 6205 in the cen- 
tral row. The crucial point is that the relative position of different 
ridge lines changes depending on the colour index adopted. This is 
clearly seen for NGC 6121, where the green line is bluer than the 
magenta line in U — B, but the opposite occurs in the other two 
planes. A similar trend is seen in NGC 6205, where the bluest se- 
quence in U — B (the green line again) is the reddest in the other 
planes. To quantify this effect, we calculate the colour difference 
between the ridge lines, at a fixed magnitude level roughly corre- 
sponding to a point two V magnitudes brighter than the MSTO. In 
particular, we fixed V = 14.82 mag for NGC 6121 and V=16.63 



mag in the case of NGC 6205. We calculated the colour differences 
with respect of the green line: the values are reported in each panel, 
and we plotted them in the bottom panels of the figure as a function 
of the central wavelength of the adopted filter. Namely, we show 
the A(t/ - S), A(S - V), and A(B - /) differences. The evidence 
that the ridge lines corresponding two the distinct sub-populations 
change their relative location in the various observational planes is 
reflected by the fact that the colour difference changes sign. From 
this plot one can immediately conclude that the colour difference 
({/ - B) — (B — I) is that one that maximizes the split between the 
different RGB sequences. Also note that the Cu,b v = (U-B)-(B-V) 
combination can work, but less efficiently than the cub.i index. We 
shown an example of this in the Appendix. 

We note that the cub.i index provides a pseudo-CMD mor- 
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Figure 13. Acu.b.i histogram for the 22 clusters analyzed. The broad or complex morphology with multiple peaks is common to all globular clusters 



phologically similar to the cy index introduced by lYong et all 
(2008), that combines three Stroemgren colours. The advantage 
of the c'u b.i index is, from the observational point of view, that 
it is based on broad-band photometry, more accessible and less 
time consuming than th e narrow-band technique. As shown by 
lYong & Grun dahl (2008), the cy index is a good tracer of the abun- 
dance of N in RGB stars. As for the case of cy, the reason for the 
efficiency of the c v B j index to split subpopulations along the RGB 
has to be searched in the spectral properties of the different popula- 
tions. This will be discussed in 5j6] However, we stress here that an 
efficient photometric identification of different stellar populations 
along the RGB is the first step to guide future spectroscopic cam- 
paigns. 



5.2 The 

c u,b,i index: multiple RGBs in 22 GCs 

The M v vs c'u b,i diagrams for a sample of 22 clusters observed in 
this survey are shown in Fig. 161 41 llR To better compare the differ- 
ent clusters, we have adopted the values of distance and average 



The cu,b,v diagram for NGC 6712 is presented in the Appendix. 



reddening from (Harris 1996, updated as in 2010) and shifted the 
diagrams to absolute magnitudes. 

Clearly, the RGBs present large colour spreads, and in most 
cases multiple sequences are evident. This strongly suggests that 
multiple stellar populations are a common feature in globular clus- 
ters. Note that the clusters are sorted for decreasing metallicity in 
Fig.loUm By comparing the pseudo-CMDs one can disclose inter- 
esting features. First of all, the slope of the RGBs in this diagram 
is strongly sensitive to metallicity. In fact, while the cu.bj index of 
the RGBs in the most metal-rich clusters ([Fe/H]> -1.3) gets less 
negative for brighter magnitudes, the more metal-poor the cluster, 
the more negative in index for the brightest stars. For metallicities 
close to [Fe/H] 1.3 the RGBs are almost vertical. 

Qualitatively, not only the slope but also the morphology of 
the RGB change from cluster to cluster. In particular, some cluster 
display a RGB dominated by two sub-branches, the best examples 
being NGC 6981 and NGC 5904, while others are characterized by 
a smoother distribution of stars (see NGC 2808, NGC 6205). Nev- 
ertheless, continuous does not imply homogeneous distribution. 
Overdensities and substructures are evident, suggesting that in all 
clusters there are dominant subpopulations, but in none of them the 
RGB appears neatly separated into several well defined sequences. 
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Figure 14. The width of the RGB, Wrgb, is shown as a function of a variety of cluster global parameter. In each panel we report the Pearson correlation 
coefficient. In case of My and iron content - shown in the top panels - we also present a linear fit to the sample (dashed gray lines) and the relative coefficients 
(al,bl,crl). The black lines show a similar fit once removed the two outliers, NGC 6366 in the case of My and NGC 7078 in the case of metallicity. 



The presence of a few stars in between might be due to photometric 
error, but we cannot exclude that they belong to tails of the distri- 
butions, with intrinsic intermediate properties. 

5.3 Dependence on the clusters'properties 

To further investigate the influence of metallicity and other cluster 
parameters on the multimodal RGBs, we have adopted the proce- 



dure illustrated in Fig.[l2]for NGC 6205. Panel (a) is a reproduc- 
tion of the Right panels of Fig. [4] showing the V versus cu,r,i di- 
agram zoomed around the RGB, where we have superimposed 
the red-dashed fiducial line, that is drawn by hand. The two gray 
lines are located 1.5 and 4.5 mag (in the V band) above the MS 
turn-off and delimit the region where the RGB multimodality is 
better visible. Only stars in this magnitude intervals are used in 
the following analysis. The normalized V versus Acu,b,i diagram 
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is plotted in panel (b), and is obtained by subtracting from the 
c u,b,i colour index of each stars the corresponding value of the 
fiducial at the same V mag. Panel (c) shows the Ac UB ,i histogram 
distribution. The histogram distributions for all the clusters studied 
in this paper are plotted in Fig.[l3]for stars in the same magnitude 
interval. Superimposed to each of them, we plot the kernel den- 
sity distribution, obtained using a gaussian kernel. The Cu,b,i dis- 
tribution of RGB stars confirms the multimodality detected in 
the pseudo-CMDs. NGC6366, NGC 104, NGC5904, NGC6121, 
NGC 6254, NGC 288, and NGC 6981 have clearly bimodal RGBs, 



while NGC 2808, NGC 6205, NGC 6218, and NGC 6254, present 
some hint of three components. 

Starting from the normalized colour histogram previously de- 
fined, we defined the width of the RGB, Wrgb, as the colour ex- 
tension of the histogram after removing the 5% bluest and red- 
dest stars (as an example, we see the two vertical lines in panel 
(c) of Fig. [120. The value of the Wrgb index, summarized in Tab. 
[2] spans a limited range, from 0.06 to 0.20 mag, with mean value 



2 The outliers have been rejected to minimize the effect of field-star con- 
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Figure 16. Acrj .b,i vj [O/Fe] abundance anticorrelation, found in all the clusters studied. 



< Wrgb >= 0. 13 ± 0.03 mag. We searched for possible correlations 
between the RGB width Wrgb and the global clusters'parameters. 
We adopted the compilation by iHarrisI (I2010T) . and calculated the 
Pearson linear coefficient with a number of different parameters. 
Some examples are shown in Fig. [14] The three rows respectively 
show, from top to bottom, the trend of Wrgb as a function of clus- 
ter integrated absolute visual magnitude M v and the iron content, 
the ellipticity and the concentration, the internal velocity dispersion 
and the central density. The Pearson coefficient is reported in each 
panel. In general, correlations seem to be very weak. The only two 
exceptions are the integrated absolute magnitude and the metallic- 
ity. Using the full set of clusters, we derive correlation coefficients 
r=-0.286 and 0.768, respectively. In both cases, plot the linear fit 
together with the data points (dashed grey lines) and report the fit 
coefficients (zero point al, slope bl, and dispersion <rl). Note that 
in both plots there are outliers with respect of the general trend. 
In the case of M v , Wrgb seems to decrease for increasing (fainter) 
magnitudes, but the faintest cluster in the sample, NGC6366, is 
an obvious outlier. If we remove it, the correlation improves from 
-0.286 to -0.619. The fit without this cluster gives a steeper rela- 
tion (black line, with a2, b2, <t2 coefficients). Note that NGC 6366 



tamination and objects with large photometric errors on the estimate of 
Wrgb. 



is a bulge cluster, and we cannot exclude that the low luminosity is 
due to severe mass loss during its orbital motion. 

A similar situation occurs for the dependence on the metallic- 
ity. While there is a general trend of decreasing Wrgb f° r decreas- 
ing metallicity, the most metal-poor cluster, NGC 7078, presents a 
colour spread larger than expected according to that relation. Re- 
moving this object, the correlation increases from 0.768 to 0.850, 
resulting again in a steeper relation (black line). A possible ex- 
planation is that the cu b.i colour becomes less sensitive to the 
colour spread of the RGB in this metal-poor regime. This could 
be due to the fact that more metal-poor RGBs are also bluer and 
hotter, and therefore the spectral features responsible for the colour 
spread might be weaker. However, it is worth notin g that a simi- 
lar ano malous behaviour of NGC 7078 was found by Pancino et alj 
feoich . when considering the CN-CH anticorrelations as a function 
of metallicity for 15 clusters. 

Interestingly, we stress that solid evidence supports the exis- 
tence of multiple populations in both NGC 6397 (the cluster with 
the smallest Wrgb , and NGC 7078, coming from both sp ectroscopic 
dLind et alj|201ll) and photometric dMilone et ai]|2012ch investiga- 
tions. In particular, NGC 6397 presents small variation of light ele- 
ments, in particular Na and O. This might also explain the limited 
colour spread of the RGB we reported here. Interestingly, when 
comparing the RGB morphology in the V vs Cu,b,i plots, one could 
speculate that the small colour spread observed in NGC 6397 is due 
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Figure 17. [Na/Fe] vs Acu b.i correlation, with overplotted the coefficients of the linear fits. 



to the absence of the red RGB, which is populated by the Na-rich, 
O-poor population. This is in nice agreement with the spectroscopic 
measurements, suggesting that the cu.b,i colour can be not only 
an effective instrument to detect multiple populations, but also to 
constrain the relative importance of different sub-populations in a 
cluster. 



6 DISCUSSION: THE CHEMICAL PROPERTIES OF 
MULTIPLE SEQUENCES 

To demonstrate the effectiveness of the cu.b.i index, in the follow- 
ing we will take advantage of the measurements of chemical abun- 
dance available in the literature for NGC6121 and NGC6205, to 
get information on the chemical composition of their multimodal 
RGBs. 

Figure Q3] is a reproduction of the right panels of Figure |4] 
zooming on the RGB of NGC6121 (top) and NGC6205 (bottom) 
in the V vs cu.b.i plane. The colour code is the same, as well as the 
definition of RGBa,b and RGBa,b,c for the two clusters. In case 
of NGC6121, s tars with avai l able s pectr oscopic measuremen ts of 
Na and O from lMarino et alj d2008l) and lCarretta et al.l d2009l) are 
plotted as circles and triangles respectively, while the Na-0 anti- 
correlation determined by these authors is shown in the two insets. 
A visual inspection immediately reveals that stars from the RGBa 



and RGBb are clustered around two distinct values of Na and O. 
The RGBa has a chemical composition similar to that of halo stars, 
and the RGBb is populated by O-poor and Na-rich stars. Similar 
conclusions can be als o inferred from the lo wer inset, showing the 
Na-0 correlation from lCarretta et al1d2009l) . even if the separation 
between RGBa and RGBb stars is less clear, likely due to the poorer 
accuracy of this data set. 

In the case of NGC6205, we identified three RGB compo- 
nent s. The Na-0 anticorr elati ons for stars of this c luster is available 
from lCohen et all J2OO5I) and lSneden e t al. (2004), and is plotted in 
the inset. Stars with both spectroscopic and photometric measure- 
ments are marked with coloured circles and triangles. We note that 
the three groups of RGB stars have, on average, different abun- 
dances of Na and O. RGBa stars are Na-poor and O-rich, RGBc 
stars are Na-enhanced and O-depleted, and RGBb have intermedi- 
ate properties. We emphasize here that the choice of the two or three 
groups of stars is not intended to demonstrate that these groups 
correspond to discrete stellar populations. The intent here is to in- 
vestigate the connection between multimodal RGB and chemical 
composition. 

Figures [T6] and [T7] display Ac v B1 versus [O/Fe] and the 
[Na/Fe] versus Acu.b.i for the seventeen clusters with available 
spectroscopic measurements. Note that for some clusters various 
spectroscopic analysis are available. In these cases we only report 
the source which provides the smallest scatter in the [Na/Fe] ver- 
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Figure 18. Panel a: Reproduction of the Actj,bj-0 anticorrelation of Fig. 16 forNGC 104. Green, mage nta, and blue circle mark the three populations a, b\ , and 
bi. Panel b: Na-O anticorrelation. Red segments separate the three groups of P, I, and E stars defined bv lCarretta et"afl f2009). Panel c\ : Na-Acu b,i correlation. 
V vs cij,b,i diagram. 



sus Ac U B i relation. Fig.[l6]shows a clear anticorrelation between 
Actj,b,i and [O/Fe]. Remarkably, a similar relation between the 
Stromgren index cy and the nitro gen abundance h as been observed 
for RGB stars of NGC 6752 bv lYong et all < l2008h . 

Interestingly, there appears to be a linear correlation 
[Na/Fe]=aAc U B i+b. The values of a and b change from one clus- 
ter to the other and are shown in Fig. [T7] together with the scat- 
ter a around the fit. Apparently, the slope of the [Na/Fe] versus 
Actj,b,i relation is different in the different clusters, and the values 
are reported in Tab. [2] 



Some hint of Acu3j-Na correlation and Acu B I -0 anticorrela- 
tion is visible also in NGC 6397, despite the small Cu.b.i broaden- 
ing observed in this cluster. This supports the presence of multiple 
stellar generations. Overall, these results imply that the c VSi index 
can be used to trace the overall distribution of sodium, and oxygen 
in GCs. 
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6.1 New hints on the stellar populations of NGC 104 

An accurate analysis of the stellar populations in each individual 
GC is beyond the scope of this paper. The study of NGC 104 (47 
Tuc) presented in this section is only meant to emphasize, again, the 
effectiveness of the Cu,bj index. As already mentioned in Sect. 12.21 
a number of recent studies have focused on the multiple stellar pop- 
ulations in this cluster, including works based on HST photometry 
and high-resolution VLT spectroscopy. In the following, we show 
that the ground-based photometry of the SUMO dataset can pro- 
vide a view on the stellar populations in NGC 104 that is even more 
complex than that given by HST. 

To do this, we start the discussion from the RGB, where HST 
and ground-based photometry revealed two distinct sequences that 



correspond to two stellar populations dMilone et al.ll2012el) . The 
less-populated component (population a) contains about 30% of 
the stars and has a chemical composition similar to the halo-field 
stars. The most-populated one (population b) includes the remain- 
ing ~70% of stars, it is enhanced in nitrogen and helium, and de- 
pleted in carbon and oxygen. 

The Acu,B,i-[0/Fe] anticorrelation of Fig. I18fa confirms the 
presence of two stellar groups. A population of O-rich stars with 
low Acu b.i values, and a population of O-poor stars with large 
Actj,b,i • These two st ellar groups correspon d to the populations 
a and b identified by iMilone et alj d2012eh . Interestingly, while 
stars belonging to population a are fairly homogeneous in [O/Fe], 
population-^ stars span a wide range of about 0.6 dex in oxy- 
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Table 2. List of the derived RGB width, Wrgb, and of the slope of the Acu.b.i vs [Na/Fe] correlation, b. 



Cluster Wrgb Acu.b.i vs [Na/Fe] 

b [mag] r 



NGC 104 


0.19 


0.48 


0.71 


NGC 288 


0,14 


0.50 


0.79 


NGC 362 


0.13 








NGC 2808 


0.20 


0.53 


0.79 


NGC 3201 


0.10 


0.34 


0.60 


NGC 4590 


0.07 


0.35 


0.38 


NGC 5904 


0.12 


0.42 


0.75 


NGC 6093 


0.14 






NGC 6121 


0.14 


0.37 


0.87 


NGC 6205 


0.12 


0.46 


0.90 


NGC 6218 


0.13 


0.46 


0.81 


NGC 6254 


0.13 


0.37 


0.73 


NGC 6366 


0.19 






NGC 6397 


0.07 


0.14 


0.40 


NGC 6541 


0.11 






NGC 6681 


0.12 






NGC 6752 


0.10 


0.56 


0.93 


NGC 6809 


0.12 


0.60 


0.67 


NGC 6934 


0.12 






NGC 6981 


0.11 






NGC 7078 


0.13 


0.51 


0.44 


NGC 7099 


0.10 


0.44 


0.37 



gen abundance. Hence we have further arbitrarily divided popu- 
lation b into two subgroups, denoted as b t and b 2 . Population a-, 
b[-, and £>2-stars are displayed in Fig. 1 1 8fa with green, magenta, 
and blue colour, resp ectively. The Na and O estimates come from 
ICarrettaetall ( l2009t) . 

These authors suggested a way to separate the cluster mul- 
tiple stellar populations on the basis of their position in the Na- 
O plane. They defined as primordial (P) component all stars with 
[Na/Fe] ratio in the range between [Na/Fe] min and [Na/Fe] min +0.3. 
Here [Na/Fe] min is the minimum value of the ratio [Na/Fe] ratio 
estimated by eye. The remaining stars are all considered second- 
population^) stars, and have been further divided into two groups. 
Stars with ratio [0/Na]> -0.9 dex belong to the intermediate (I) 
population, while those with [0/Na]< -0.9 dex are defined as ex- 
treme (E) population. E, I, and P groups are separated by red lines 
in the [Na/Fe] vs [O/Fe] plane shown in Fig.|18b. The P component 
includes all the stars in population a, but is strongly contaminated 
by population by stars. Similarly, the I component contains a mix 
of population bi and population b 2 stars. 

An effective way of testing the presence of a real oxygen 
spread amongst population b stars, makes use of a comparison 
of independent measurements of different quantities for the same 
stars, like Ac UB I , [Na/Fe], and [O/Fe]. Specifically, if the broad- 
ening is intrinsic, population a-, b\-, and £>2-stars, will also differ 
in their sodium abundance. But if the broadening is entirely due to 
observational errors, the stars in the three stellar populations have 
the same probability of being Na-rich or Na-poor. Panel Ci of the 
same figure shows Actj.b.i against [Na/Fe], and the histogram dis- 
tribution of [Na/Fe] is plotted in panel c%. The fact that popula- 
tion a-, b]-, and fc-stars have, on average, different sodium content, 
demonstrates that the oxygen spread is intrinsic and that population 
b contains further sub-populations. The position of stars belonging 
to the three stellar populations in the V versus cu.b.i diagram is 
shown in panel (d). 

Independent support to the fact that the star formation his- 



tory of NGC 104 is more complex than what inferred from pre- 
vious works comes from the analysis of HB stars. Both theo- 
retical arguments and observations have shown that, at least in 
some clusters, the HB morphology is st rictly related to the pres- 
ence of multiple s t ellar populations (e . g. |p ' Antona & Calo feooljl ; 
ISalaris et alJl2008l : lMarino et al.ll201 lUGratton et al1l2012l) . In case 
of NG C 104, there are two gro ups of CN-rich and CN-poor red-HB 
stars dNorris & Freeman! 1982b . that define two sequences along the 
red HB, when observed in appropriate filters (Milone et al. 2012). 

Figure [T9] shows data for NGC 104. We plot in panel (a) a 
CMD based on wide colour baseline, with the sample of HB stars 
shown as black dots. In panels (bl) and (b2) we show the V ver- 
sus Ctj b i CMD and Hess diagram, respectively. The correspond- 
ing histogram of the colour distribution is shown in Panel (b3) and 
confirms that the HB of NGC 104 is not consistent with a simple 
stellar population. There are at least two groups of stars in close 
analogy with what observed along the RGB. A clustered blue HB, 
and a more dispersed red one with the presence of at least two com- 
ponents. Interestingly, stars in the red component are also brighter 
than the stars in the blue portion of the HB. Moreover, a similar 
conclusion can be reached also for another cluster with predom- 
inantly red HB morphology, namely NGC 6366 (see Fig. |6), that 
displays a bimodal HB made with similar morphology. 

Though we cannot exclude that mass loss play a subtle role in 
shaping the complex HB morphology, it is difficult to understand 
why this effect would manifest in certain colours rather than in all 
CMDs. Nevertheless, we conclude that these results on NGC 104 
provide strong direct evidence that the red-HB morphology of this 
GC is strictly related to hosted multiple generations, and suggest 
that the multiple sequences discovered in the CMDs of other GCs 
studied in this paper may be related to the morphology of their HBs. 

Figure 1 1 9b also reveals a large number of AGB stars, are 
marked with black asterisks. To investigate the presence of mul- 
tiple stellar populations along the AGB of 47 Tuc, panel (cl) shows 
the V versus Cu,b,i diagram zoomed around the AGB, with a fidu- 
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cial line through the redder component. We verticalize the AGB 
stars in the usual way in panel (c2), while panel (c3) shows the his- 
togram of the Acu.b.i distribution. The bimodal distribution demon- 
strates that also the AGB of 47 Tuc hosts at least two stellar pop- 
ulations. We emphasize that this is the first detection of multiple 
sequences along a cluster AGB. Noteworthy, NGC 104 is the only 
cluster where multiple sequences have been found from the MS, to 
the RGB, HB, and finally alo ng the AGB. 

It has been suggested by Norris (1981) that there should exist 
a tight correlation between the incidence of AGB stars belonging 
to helium-enhanced stellar generations and the HB morphology. In 
particular, the observed paucity (if any) in GCs of CN-strong AGB 
stars could be related to the fact that these objects would be the 
progeny of He-enhanced HB stars - those stars populating the ex- 
tended blue tail of the HB in some GCs - and given these stars 
would evolve preferentially as AGB-manque objects as a conse- 
quence of their extremely thin envelopes, they would miss the AGB 
stage. In this context, the ability of c U B j to trace the presence of 
different sub-populations also along the AGB appears useful in or- 
der to investigate this issue in a large sample of clusters. 



of view of the ground-based photometry, will allow us to investi- 
gate fundamental aspects in our understanding of multiple stellar 
populations, e.g. the radial distribution of different stellar groups 
and the real fractions of stars belonging to different stellar com- 
ponents. Finally, we stress that other color combinations can be 
used to distinguish multiple populations in the CMDs, like c V B y — 
(U — B) — (B — V). However, the separation between the different 
sequences is smaller than that provided by cu b.i index. 
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7 CONCLUSIONS 

This is the first of a series of papers that will present the results of 
the SUMO project - a SUrvey of Multiple populations in globular 
clusters. We focused here on a general description of the project, 
and the first general results. We are building a large data base of 
homogeneous photometry of a large number of GCs, with the main 
aim of characterizing and interpreting their multiple populations. 
This will be basis for future spectroscopic follow-up. To date, we 
have collected (/-band images for a sample of ~30 clusters com- 
plementing the archive maintained by P. B. Stetson. In this paper, 
we have summarized the observational campaigns realized so far, 
the data reduction strategy, and the techniques to minimize differ- 
ential reddening and polluting sources like foreground stars and 
background galaxies. 

We have introduced a new index,cy g j , defined as the differ- 
ence (U — B) — (B — I), that turns out to be very efficient in dis- 
entangling multiple sequences along the RGB in the range of ages 
and metallicities covered by selected objects. We presented results 
for 23 clusters, showing that showing that: i) they all display clear 
evidence of broadened or multiple RGBs; ii) the slope of the RGB 
in the V-Cu B j diagrams is very sensitive to the mean cluster metal- 
licity. 

Although in those clusters where multimodal RGBs are evi- 
dent, small intermediate populations always appear. These could be 
outliers due to photometric errors, but we cannot exclude popula- 
tions with intrinsic intermediate properties. This continuous but in- 
homogeneous distribution seems to disfavor formation mechanisms 
that predict bursty formation of subpopulations. 

The properties of the Cu,b,i index were exploited to revise the 
behaviour of stellar populations in NGC 104, showing an overall 
more complex picture than previously discovered using HST data. 

The effectiveness of the Cu,b,i index was validated relating 
the photometric and spectroscopic properties of RGB stars in sev- 
enteen clusters. In particular, Na-rich and O-poor stars populate a 
RGB with larger (less negative) c v B l index than the one defined by 
the Na-poor, O-rich sub-population. These results demonstrate that 
the our photometric approach is an alternative to the spectroscopic 
one in the identification of GC stellar populations. In addition to 
this, the study of all the RGB stars provided from the wide field 
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APPENDIX A: THEc U B V AS AN ALTERNATIVE INDEX 
TO IDENTIFY MULTIPLE RGBS 

Fig. IA II shows the case of NGC6712 (panels al-a4). Since the I- 
band data are not yet available for this cluster, we used the Cu,b,v = 
(U — B) — (B — V) index. Note that the sensitivity of this index is 
smaller than that of the Cu,b,i, as demonstrated by Fig. [5] In fact, 
the colour difference between difference sequences is smaller in 
(B - V) than in (B - I). Nevertheless, it allows to identify well- 
defined substructures with at least two RGBs. To support this con- 
clusion, we compare the Cu,b,v index of NGC 6121 (panels bl-b4), 
which presents clear hints on multiple RGBs from the c VBl index. 
In particular, panels al and bl present a zoom in the RGB region 
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with superimposed the ridge line of the reddest sequence, in anal- 
ogy to what presented in Fig. [T2] The rectified RGBs are shown 
in a3) and b3), and finally the histograms showing the distribu- 
tion of Acubv colour differences are plotted in a4) and M). The 
latter clearly shows bimodal distribution, thus supporting that also 
NGC 6712 hosts multiple populations. 
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Figure Al. My vs cv,B,\ of the cluster NGC 6712, for which / photometry is not yet available, in comparison with that of NGC 6121. The lower panels show 
the rectified RGBs, and the Acu b.v colour distribution, which clearly show hints of bimodality. 



